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I. INTRODUCTION 



Since 1988, a non-extensive formalism of statistical mechanics jlj has been developed and up to recent 
years, many works have been devoted to show the robustness and usefulness of this approach. Wc believe it is 
robust in the sense it allows generalizations of a variety of fundamental concepts of statistical thermodynamics 
0, such that it avoid to enter in severe contradiction with well established facts. In addition, we believe 
it is useful in the sense it provides a theoretical basis and relevant explanation of some experimental and 
observational situations ||, where Boltzmann-Gibbs statistics fails. That is to say, where Boltzmann-Gibbs 
statistical functions diverge and do not yield any physical prediction. 

However, this formalism, unlike Boltzmann-Gibbs', has a non-zero set of free parameters, here represented 
by Q {q £ 3?)- This unique parameter controls the degree of the nonextensivity of the system in consideration. 
(At this point, it is worth noting that the formalism includes the standard, extensive, statistics as a special 
case for the value of q = 1 and all expressions derived within this non-extensive framework give the results 
of Boltzmann-Gibbs statistics in the q — * 1 limit.) Only after 1995, some works started to address the long 
standing puzzle of understanding the physical meaning of q. Amongst the works related to this topic, two 
main streams started to become more apparent. On one side, there are attempts on the study of conservative 
j| and dynamical systems, more precisely, dissipative systems with both, low || and high || dimensions. 
On the other, there have been some efforts of estimating bounds upon q in measurable physical systems. The 
physical applications studied so far can be enumerated as follows: the microwave background radiation @J^], 
the Stefan-Boltzmann constant |^,|l(| , the early Universe |ll],[l2| and the primordial neutron to baryon ratio 
in a cosmological expanding background |L3]]. As a final note, a recent letter by Alemany jl4|] explored the 
definition of a new fractal canonical ensemble, associated with the parameter q. However, first estimations 
for the universe as a whole yielded to values of q bigger than those allowed by nucleosynthesis |b|] . 

In all these works, two different approaches for the quantal distribution functions have been used: a closed 
analytical form for them is still lacking. The asymptotic approach of the kind (3(1 — q) — > of Tsallis et 
al. fl7| was used in Refs. @,|,|llp| and also in some other applications such as the study of Bose-Einstein 
condensation in a fractal space [ lM , the specific heat of 4 He |l6| and the thermalization of an electron-phonon 
syste m fL7[ . The second approach for the generalized distribution functions has been proposed by Buyiikkihc 
et al. pW, which we refer to as Factorization Approach. This last term is justified in that it is generically 
based on the factorization of the generalized grand canonical distribution and the concomitant generalized 
partition function as if they were extensive quantities. However, although the generalized distribution 
functions within factorization approach have been derived before the asymptotic approach, they have not 
been prefered to be used in physical applications, mainly due to a work by Pennini et al. |p_9|| . There, the 
authors claimed that the quality of the results of the factorization approach deteriorates as the number of 
the particles of the system increases. But very recently, contrary to the previous belief, new results by Wang 
and Le Mehaute [ p0[ , favoured the factorization approach and showed clearly that there exists a temperature 
interval where the ignorance of the approximation is significant, but otherwise, the results of this approach 
can be used with confidence no matter the number of particles. In fact, they showed that, for a macroscopic 
system (i.e., with ~ 10 23 particles) having two states with a small energy interval of about 10 10 Hz, this 
forbidden temperature zone is very narrow and situated at extremely high temperatures (~ 10 20 K) and as 
the number of levels increases the forbidden temperature zone shifts to higher temperatures. Therefore the 
results of the factorization approach can be used at low temperatures up to 10 20 K without any ignorance 
for any physical system under consideration pc|-|22]| . As it was expected, the study of Wang and Le Mehaute 
started to accelerate new attempts of ap plying these distribution functions to physical systems in order to 
estimate some alternative bounds upon q |lO|Jl^ ] ; being this approach more handable than the one advanced 
by Tsallis et al. Although, after Wang and Le Mehaute |^0), some physical systems have been worked out 
with the help of the generalized distribution functions of the factorization approach, general expressions for 
magnitudes of fermions and bosons are still lacking. 

Indeed, the possible need for a nonextensive formalism of thermostatistics is clear for a long time in 
gravitation 23 1, magnetic systems [Q, Levy- like diffusions |2f|, some surface-tension problems [ p6[ , etc., 
since Boltzmann-Gibbs formalism is known to fail whenever the physical system under consideration includes 
(i) long-range interactions and/or (ii) long- memory effects and/or (iii) the system evolves in a multifractal- 
like space-time. The last property can be understood as the system evolves in a porous-like medium where 
the properties of the multifractal structures govern it. 

In the light of the facts stated above, we address the study of the generalized quantal distribution func- 
tions of fermions and bosons in a non-Euclidean, fractal-like space-time and to obtain general results for 
magnitudes, such as the number of particles and the internal energy that could be used in general for 
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any application to generalized systems. Once these general formulae be established, we proceed further 
investigating how some simple constructs of the statistical theory are affected by the change of distribution 
functions, for instance, the Chandrasekhar mass and the Bose-Einstein condensation. A similar work study- 
ing the Bose-Einstein condensation was recently done, as commented above, by Curilef |l5|] and our results 
are compatible in this case. Finally, it is worth noting that, for non-interacting fermions and bosons, the 
information about the fractal structure is kept in the nonextensivity parameter q, which is by now verified 
to be connected with the fractal dimension and the multifractal singularity spectrum ||. In a sense, we 
shall study if any correction to the usual number of particles and energy arises due to the change of the 
distribution functions. 



II. THE FERMION CASE 



Let us start by writing down the generalized distribution function for fermions, up to (1 — q) order. It is 
given by |J), 



hlf- 



ermions\ 



1 q - 1 [f3(e - /i)] 2 e^ 6 -^ 

2 /3(e- M ) + 1 + ~ ( e /3(e-M) + 1) 2 ' 



(1) 



where (3 = 1/kT, e is the energy level and fx is the chemical potential. 

We are interested in solving the following integrals. For the number of particles, 



<N>= 



4irV(2m) 3 / 2 



de 



d/2 



9-1 



e /3(e-p) + I 



00 d (/?(e- ^fe^-^e 1 ' 2 



(e0(*-f) + 1)2 



(2) 



and a similar one, replacing the factors e 1 / 2 for e 3 / 2 , for the energy < U >. Note that we have assumed the 
same pre-factor for both terms. This can be justified if one thinks of this factor as coming from microscopic 
quantum considerations and a change of variables, from momentum to energy, inside the integrals. See for 
instance, section 2.4 of Ref. p7[ . 

Then, we can split < N > into two parts. Firstly we have the usual expression, given by the first term in 
the rhs of (S), 



U sual = 



h 3 



8TrV(2mfi) 3 / 2 



fite-n) + 1 



3h 3 



1 TT 



7 f TT V 



640 V/W 



And we also have a second term, 



_ 47T^(2m) 3 / 2 q - 1 f°° ((3{e - ^fe^^e 1 ' 2 _ 4ttT/(2to) 3 / 2 q-l 
q ~ P 2~ J € + 1)2 ~ h 3 2~ ' 

I can be recasted using the following dimcnsionless variables: 

x — /3e, ip = Pfj,. 

This enables us to write 

- ijjfe^x 1 / 2 



I = P 



dx = (3- 3/2 I x . 



(e x -^ + l) 2 

In order to solve this last I x we apply the following trick. We know that integrals of the form 

m 



dx 



(3) 



(4) 



(5) 



(6) 



(7) 



have a solution, 
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W , ^ df , 77T 4 d 3 / 317T 6 d 5 / 

J W + 6 dx 360 dx 3 15120 dx 5 



(8) 



where the derivatives have to be evaluated in x = ip p7[ . In fact, this was already used to write down the 
solution for the usual term of (||). The idea is then to introduce an extra parameter, Z, in the integral ([?]) 
and consider, 



fix) 



dx. 



(9) 



Here, it is important to stress that f(x) does not depend on I. If we derive inside this integral with respect 
to I, the result will be 



f(x)ipe 



x—lifj 



rdx, 



lo (e^' + l) 2 ' 

and if / = 1, it can reproduce our integral I x iff the function f(x) is given by, 



/(*) 



Thus, using (O), the solution to I x is, 



L = 



(X - VP 2 t 

— ; — x 

' fix) 

dl J (e x -^ + 1) 



/2 



dx 



(10) 



(11) 



(12) 



i=i 



We have to be very careful in using the correct formula to get the integral in fix): we should redefine ip = lip 
and apply then (ph. So, the mechanism is, considering f(x) as in (O), compute the integral (||), derive the 
result with respect to I and finally evaluate in I = 1. That yields, as we have seen, the integral we need.[] 
Let treat each of the term in m) separately. The first term is obtained inmediately: 



which finally leads to, 



jx - iP) 2 ! 

— / — x 

Ip 



1 st term 



1 f$ 

/2 dx = - / (x 2 - 2xip + ip 2 )x 1/2 dx, 
W Jo 



1 

7' 



2lp^— 



7/2 ' 5/2 v 3/2 
Now we have to derive with respect to I (recall that ip = lip), obtaining 



(14) 



(15) 



d 1 st term 1 
dl ~~ ip 



( 7 /2)l 5 / 2 ^— - 2(5/2)? 3 / 2 f TT + (3/2)Z 1 / 2 ^ 



and evaluate at I — 1. The result is, 



(7/2) 



d I s * ierm 
dZ 



(5/2) 



(3/2) 



(16) 



(17) 



1=1 



1 A more transparent way of solving these integrals (see acknowledgments) is to take the integral (6) as 

■st / ; — rdx-2ip— / , — -di + f — / , — -dx (13) 

dip J Q e x -^ + l Y dip J Q e x -^ + l Y dip J Q e x -^ + l y ' 

and apply expansion (8) to each part. As we see, in fact there is no need to introduce such an extra parameter I, 
and it only stands as a mathematical trick in a particular form of solving the integrals. Numerical results for the 
coefficients obtained with both methods are the same. We warn the reader not to be distracted for such a trick. 
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The second term in M) is obtained as, 



2 nd term = — 
6 



i> 2 



which finally yields to, 

2 nd term = 

The derivative with respect to I is, 
6 d2 nd term 



7T 



2 vr ■ ry W 



Evaluating in / = 1 we obtain, 



+ ^-^^-3/2^-3/2. 



d2 nd term 



1/2 



i=l 



One can do the same with the third term. In this case we need to compute the third derivate. It is 

-3/2 



r 



o -r-3/2 o ™-5/2 

^- 1/2 -^-«V +(i -«lV- 



Computing this last for x = ip and deriving with respect to I we get, 



dl 



-15 
To" 



(ty) 



-7/2 



And finally making 1 = 1, the correction is found to be, 



dl 



J 2 = 1 



Now we can go back and collect the results for /, which results, 

/ = p^I x = 0-3/2 ( + — - —i>- Zl2 

' H \ 3 120 



and the final expression for < N > is then, 

3h 3 



<N> 



1 fit 



8irV(2mfxf/ 2 



7 / rr 



8 V/W 640 V/3m 



9-1 







rr 2 1 21rr 4 1 



2 (3li 240 (/3^) 3 



This final expression represents the correction terms due to non-extensivity that arise for < N >. 
We can do exactly the same for the energy < U >. I q is given now by, 



_ 47ry(2m) 3 / 2 g-l 5/2 
9 ~ ft3 2 ' 



where 7 X is, 



2* = 



(a - VQV^^ 

(e x -^ + l) 2 



We apply again the same trick, in this case, with a function g(x) — (1/i/j)(x — ip) 2 x 3 / 2 . As before, 



9{x) 



(e x -W + 1) 



dx 



(29) 



1=1 



The first term in the serie (|J), is again equal to zero, as can be directly verified by computing the integral 
for g, deriving the result with respect to I and evaluating at I = 1. The second term needs 



dx 



1 



Derivating this with respect to I and evaluating at I = 1 we obtain for the correction, 

^2 



3/2 



The third term needs the third derivative of g. This is found to be, 



1 



Ib- 



-1/2 



{x - ij))- - {x - tpf -x 



-3/2 



Again deriving with respect to I and evaluating at I = 1, it yields the following correction 

7^ 4 



360 



9V> 



-1/2 



(30) 



(31) 



(32) 



(33) 



Collecting the results, we get for < U >, 

5h 3 

< U > 



87ry(2m) 3 /V/ 2 



- 1 



-V 

8 V/W 
5tt 2 1 

0+ — w 



7 / 7T 

" 384 [Jjj, 

315tt 4 1 
720 (/3 M ) 3 



(34) 



As above, this expression represents the correction due to non-extensivity that arise for the energy. 

It is important to note that the T — > limit is attained without any correction thus suggesting that 
non-extensivity can hardly have any role for very low temperatures. One can now study particular physical 
systems. For instance, one can inmediately realize that the Chandrasekhar mass will be not modified at all 
in passing to a non-extensive context, at least up to order (1 — q). This is simple because the Chandresekhar 
mass is a construct that arises at T = where both corrections, to < N > and to < U >, are exactly 
zero. However, if non-extensive statistics is concerned in the study of stellar structures, corrections will 
unavoidably arise. 



III. THE BOSON CASE 



Now we analize the boson generalized distribution function, namely, 

1 q- 1 [f3{e - (i)] 2 e^-v) 



^q[bosons] 



( e /3(e-M) - 1)' 



(35) 



where again (3 = 1/kT, e is the energy level, fi is the chemical potential. Then, using equation (|35|), one can 
obtain the average number of particles as 



(N) 27r(2m) 3 / 2 



V 



°° e 1 ' 2 de q-l f°° /3 2 (e - ^) 2 e^-^e 1 ' 2 de 



o (e^(«-/*) - 1)* 



(36) 



As in the standard case, we have to separate the state with e = 0, which has zero weight in the integral (pq). 
For this level of energy, the distribution function yields, 
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z q — 1 flnz) 2 ! 

^ o )~[ i+ Vt-h ^ 

where z = is the fugacity of the gas. In the case z <C 1, the correction term goes to zero as z(lnz) 2 does. 
If z —* 1, the correction goes to 1 as z, and results are neglectable in comparison with the first, diverging, 
term. This can be seen in Fig. 1. 

Substitution of x = f3e and ip — /3(i in equation (Bq) yields, 



(N e ) q 2n{2mkfl 2 



rp3/2 



(38) 



V h 3 

where (N e ) stands for the number of particles in the excited states (e ^ 0), and 1\ and I q are defined to be 

f°° x 1 ' 2 dx f°° (x - ^) 2 e x ^'x 1 / 2 dx 

h -Jo ^-J (e ,-* _ l) 2 ■ (39) 

The integral 7i is the one which appears in the standard, q — 1 case of the boson gas and therefore the 
solution of it is the known one p7[ , 

h =r(3/2).g 3/2 (2), (40) 

where <7n(z) = y^y /o°° z-^e^-i ^ 3 ' — EJli ( z ' s /s n ) f° r small z and T is the usual Gamma function, T(n) = 
Jo" s~ x x n ~ 1 dx. For < z < 1 and Vn, n > 1, the functions g n {z) are bounded by the Ricmann zeta 
functions, which yields for all z values of interest, 

h < r(3/2)C(3/2). (41) 

On the other hand, the integral I q is again somewhat cumbersome to solve, but after some algebra, similar 
to what we did in the previous section, we managed to find the analytical solution as 

I g = T(7/2)g 5/2 (z) - 2VT(5/2)«? 3/2 (z) + T(3/2)^ 2 g 1/2 (z). (42) 

To get the previous result one has to take into account the known relationship between the g n (z) and its 
derivatives, 

d d 

g n -i(z) = z—g n (z) = — — -g n {z). (43) 
oz a(lnz) 

Recalling again the definitions ip and z, we can write i\> = lnz. Putting these solutions of the integrals in 
equation (|38|), one can obtain 

(N e ) q 2*(2mky/ 2 3/2 

r(3/2). 93/2 (z) + i^l [r(7/2) 55/2 (z) - 2r(5/2) lnz 53 /2(^) + r(3/2)(lnz) 2 3l/2 (z)] | , (44) 

which is the g-dependent solution for the number of particles in the excited states of boson systems, including 
the standard case as a special one if q = 1. 

Like in the standard case, let us study each term of the correction in a separate fashion in order to find 
the bounded form of them. The first term is similar to l\ and hence it is bounded by r(7/2)£(5/2). The 
second and the third terms are different from the first one in the sense that they do not take their largest 
values at z = 1, but instead, both of them tend to when z — > and z — > 1. This behavior can be seen 
in Fig. 2. The largest values are situated at z max \ — 0.447 for the second term and z maX 2 ~ 0.175 for the 
third. 

Consequently, if all these bounds are put together, then the total number of particles in all excited states 
is also bounded by, 

< v M^f /2 Tm (r(3/2) C (3/2) + 



h 3 

[r(7/2)C(5/2) - 2r(5/2) ln(0.447) 53 / 2 (0.447) + r(3/2)(ln(0.175)) 2 9l/2 (0.175)] } (45) 
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which gives 



(N e ) q < V — — T 3/2 {2.315 +(q- 1)3.079} . (46) 



If we now concentrate on Bose-Einstein condensation, then the condition for the appearance of it can be 
expressed as 

N > (N e ) q . (47) 
Alternatively, with constant N and V, this condition can be recasted in the form, 

h 2 ( N \ 2/3 

T < {Tc)q = (2tt) 2/3 2mfc I V [2.315 + (q 1)3.079] J ' (48) 



N 



or up to order 1 in (q — 1), 

where (T c ) is the g-dependent characteristic temperature of the Bose-Einstein condensation. It is easily 
seen that this result shows that the critical temperature decreases when q < 1, which is consistent with the 
previous result of Curilef Jl5|. Note that the standard (T c ) 1 case can easily be obtained for q = 1 value in 
the above expression. 

Any accurate simultaneous determination of N, V and T can yield, in principle, a bound upon q. However 
in practice, this could be well below any practical possibility, due to the smallness of the correction term. 



IV. CONCLUSIONS 



Our main results can be summarized as follows. We have been able to solve, using the generalized quantal 
distribution functions in the factorization approach, the values for the average number of particles and energy 
in the case of system of non-interacting particles, either fermions or bosons. These results are expected to 
be useful in any analysis of statistical phenomena within the context of non-extensive scenarios. We could 
explicitly see that all the terms coming from the non-extensive part of the integrals go to zero when the 
temperature goes to zero, something that could be expected due to the form of n q . Non-extensivity can 
not play a role for very low temperatures, at least up to order (1 — q). As a consequence, for instance, 
the well known Chandrasekhar limit for white dwarfs stars is not affected by a change of the statistical 
framework. However, it is to be explicitly stated that any other model of star, that happens with T / 0, 
will be affected by such change. This is why, for example, this statistical scenario could be useful to tackle 
the solar neutrino problem |28| . In the boson case, we have seen that a small correction appears to the 
Bose-Einstein condensation, this happening not exactly at T = 0. If this correction is enough to work out 
a possible bound upon q in this system remains to be studied. In the paper by Curilef it was argued that 
due to a possible fractility of the universe, distribution functions could be modified in the sense described 
here, and that this could be useful to study the behavior of diffuse gas clouds. This could be an example of 
where the q ^ 1 statistics may arise, although due to the form in which equilibrium distributions of boson 
stars arise it is unlikely that this could be applied for those cases. 
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Figure Captions 



Figure 1 : The plot of N q (e = 0) as a function of z for various q values. 



Figure 2 : The plot of the second and third terms of the g-dependent part of eq.(43) as a function of z. 
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